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Abstract 

Midrapidity nucleoli elliptic flow is studied within the Boltzmann-equation 
simulations of symmetric heavy- ion collisions. The simulations follow a lattice 
Hamiltonian extended to relativistic transport. It is demonstrated that in the 
peripheral heavy-ion collisions the high-momentum elliptic flow is strongly 
sensitive to the momentum dependence of mean field at supranormal den- 
sities. The high transverse-momentum particles are directly and exclusively 
emitted from the high-density zone in the collisions, while remaining par- 
ticles primarily continue along the beam axis. The elliptic flow was mea- 
sured by the KaoS Collaboration as a function of the transverse momentum 
at a number of impact parameters in Bi + Bi collisions at 400, 700, and 
1000 MeV/nucleon. The observed elliptic anisotropies in peripheral collisions, 
which quickly rise with momentum, can only be explained in simulations when 
assuming a strong momentum dependence of nucleonic mean field. This mo- 
mentum dependence must strengthen with the rise of density above normal. 
The mean-field parametrizations, which describe the data in simulations with 
various success, are confronted with mean fields from microscopic nuclear- 
matter calculations. Two of the microscopic potentials in the comparisons 
have unacceptably weak momentum-dependencies at supranormal densities. 
The optical potentials from the Dirac-Brueckner-Hartree-Fock calculations, 
on the other hand, together with the UV14 + TNI potential from varia- 
tional calculations, agree rather well within the region of sensitivity with the 
parametrized potentials that best describe the data. 
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I. INTRODUCTION 



Central collisions of heavy nuclei give, principally, a unique opportunity to access the 
properties of nuclear matter at high baryon densities and at high excitation energies. In prac- 
tice, though, these collisions often represent a complex puzzle where many effects, related to 
different matter properties, compete in generating observables. If one takes e.g. the collective 
flow in central collisions, it is generally known to be affected by the dependence of nuclear 
mean field (MF) on baryon density and on particle momentum, as well as by the magnitude 
of in-medium nucleon cross-sections. Different sets of assumptions on the dependencies and 
on the cross sections could be used to theoretically explain individual data on flow. 

For the sake of progress, under the circumstances indicated above, it becomes of ut- 
most importance to isolate either regions of data or data combinations that are reasonably 
sensitive to individual rather than many uncertain features of nuclear matter, simplifying 
the reaction puzzle. Here, we concentrate on the momentum dependence of the nucleonic 
optical potential. The momentum dependence of that potential has been investigated at 
normal and subnormal densities utilizing nucleon-nucleus elastic scattering data and 
different microscopic predictions have been made [0-0] with regard to the densities above 
normal. Though supranormal densities are reached in central heavy-ion reactions, only cir- 
cumstantial evidence for the momentum dependence of the potential was seen, primarily in 



flow ||10 Hl^1 , and any quantitative assessment of that dependence was out of question. Here, 
we show that focussing on particular features of the collective flow allows not only to see the 
momentum dependence directly but, in fact, to assess it quantitatively at the supranormal 
densities. 

The momentum dependence posed difficulties for low-energy central-reaction simulations 
due to problems with energy conservation. We elaborate on recent advances that expanded 
up and down the beam energy range where simulations can be done reliably. 

In the next section, we discuss at some length the framework of our investigations. In 
Sec. 3 we draw conclusions on the optical potential by comparing results of the reaction 
simulations with specific data on the elliptic flow. In Sec. 4 we confront the findings with 
the microscopic calculations of the optical potential. We summarize our results in Sec. 5. 



II. TRANSPORT THEORY 

A. The Boltzmann Equation 

We follow the dynamics of central collisions of heavy nuclei within the Landau quasipar- 
ticle theory of which the relativistic formulation has been given in 0]. In the quasiparticle 
approximation, the state of a system is completely specified when the phase-space distribu- 
tions fx = fx{p, r ,t) for all particles are given. The distributions satisfy the Boltzmann 
equation 

dfx . de x dfx de x dfx r < n ^ , v r> , m 

The l.h.s. accounts for the motion of particles in the MF, while the r.h.s. accounts for 
collisions. The single-particle energies are derivatives of the net energy E of a system, with 
respect to the particle number: 
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, (2tt) 3 8E 
9x 5f(p,T,t) 

where gx is spin degeneracy. Vector vx = dex/dp on the l.h.s. of (|l|) is velocity, while 
—dex/dr is force. The degrees of freedom in the description are nucleons, pions, A and N* 
resonances, and, optionally, light (A < 3) clusters. Factors /C < and /C > on the r.h.s. of ([I]) 
are the feeding and removal rates. 

The net energy of a system forms, obviously, a four-vector with the total momentum. 
Since the single-particle momentum p may be represented as a derivative of the total mo- 
mentum with respect to /, similarly to (0), and / is a Lorentz scalar, it follows that (e(p), p) 
form a four-vector. This is independent of the specific dependence of e on p in any one frame. 
By considering a change of the frame, one can subsequently see that de/dp transforms as 
a velocity in the Lorentz sense, so that (7, 7 v) = u, with 7 = 1/ \/l — v 2 , form a four-vector. 

The combination of relativity and momentum dependence brings in some pecularities in 
the collision rates and cross sections, beyond what is encountered nonrelativistically. Thus, 
consistently with the Fermi Golden Rule and the requirements of covariance, the contribution 
of binary collisions of particles X to the removal rate in (|I]) is 

9x f dp 2 f dp[ f dpi 



K5(Pi) = ?T / t£^" T^tj 753b a* 



12 *-TTj—_ —12 



l2X->2X> 



7i J W12J (2tt)3 7 ^ 2' 

x (2vr) 3 5( Pl + p 2 - pi - p' 2 ) 2 vr 5(6, + e 2 - ei - e' 2 ) f 2 (1 - f[) (1 - f 2 ) 



In the above, \Ai\ 2 represents a squared invariant matrix element for scattering, averaged 
over initial and summed over final spin directions. The factors 7 are associated with the 
respective velocities and dp/7 is the invariant measure as may be verified by considering 
the change of a reference frame. The starred quantities in (^) refer to the two-particle cm. 
defined by the vanishing of the three-momentum, P = 0, where P = pi + p 2 - The cross 
section in @ is given by 

do- p*' 2 

TjtW = a <? * — *' — \M2X^2X'\ 2 ■ (4) 

dtt* Ati 2 7* 7* v* l2 7* 7* v* l2 
The relativistic relative velocity in (^) and (^) is defined through 

[{Pu^u.-iPu^u^ 112 



7l 72 Vl 2 



(5) 



p 2 

The above definitions ensure the standard form of the detailed balance relation, i.e. here 

In the cm., the relative velocity reduces to the velocity difference. The factor of 1/2 in front 
of the angular integrations in (|3]) accounts, in the standard manner, for the double-counting 
of the final states in scattering for like particles. 
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B. The Energy Functional 



The MF dynamics follows from the dependence of total energy on the phase-space dis- 
tributions. We adopt simple parametrizations for the net energy, that permit the transport 
calculations to be carried through, and that are flexible enough with regard to the MF and 
to the equation of state (EOS). In our parametrizations, the energy consists of the covariant 
volume term and of the noncovariant gradient-correction, isospin interaction, and Coulomb 
terms defined in the system frame: 

E= I dre + E 1 + E T + E coul . (7) 



In the practice of central reactions, the importance of covariance and the indispensability of 
the correction terms are usually mutually exclusive. For different systems, there is a level of 
cancellation between the Coulomb and isospin terms. The Coulomb energy in ([^D is simply 

E wal = -^Sdvdv' pA{v)pc f' ) . (8) 
4vre J |r - r'| y ' 

The gradient term is 

where p denotes the baryon density and po = 0.160 fm~ 3 is the normal density. Finally, the 
isospin interaction term is 



T 



(10) 



where px represents the density of the third component of isospin. 

The gradient term in the energy allows us, primarily, to account for the effect of the finite 
range of nuclear forces, augmanted by the lowest-order quantal effect of the curvature in the 



wavef unctions, in the Thomas-Fermi (TF) initialization of nuclei [p~6| , p~T[| for our reaction 
simulations. We take the coefficient in equal to ai = 21.4 MeV fm 2 for momentum- 
independent MFs, as corresponding to the finite-range correction from the Skyrme effec- 
tive interaction and an addition from the Weizsacker kinetic-energy term [pl|. For the 
momentum-dependent fields, we take a bit lower a% = 18.2 MeV fm 2 from our own adjust- 
ments to ground-state data. The initialization is described in the next subsection. 

The isospin term fllUD in the energy contributes to the isospin asymmetry coefficient 
in the Weizsacker mass formula (at a 50% level) and gives rise to the isospin term in the 
optical potential, of the magnitude \ N ~ Z \ for nucleon scattering off a nucleus. We use 
ax = 97 MeV, that simultaneously produces reasonable results for the mass asymmetry and 
the asymmetry in the potential deduced from nucleon scattering []l9",2fl. 



Aside from the above isospin component, the strong- interaction field (from the covariant 
volume part) is chosen as acting only on baryons in our calculations. That is the field that 
we are interested in here. Pions are, anyway, infrequent in the energy range within which 
we will be making comparisons to data. We should note that, when the vector and scalar 
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type MFs may be momentum dependent with no exclusive dependence on the vector and 
scalar densities, there is neither a benefit nor a phenomenological basis, in the absence of 
spin dynamics, for a separate consideration of these fields. 

Guided solely by the calculational convenience, we choose the fields that could be easily 



identified as either vector or scalar 10 . Thus, in the case of the fields without momentum 



dependence in their nonrelativistic reduction, we use the energy density in the form JT5[ 



52 9x 



x 



J fx(p) V? 2 + m x (Ps) + J" dpi U{p' s ) - Ps U( Ps ), (11) 



where mx(p s ) = mx + Ax U(p s ), Ax is baryon number, and 

\^ A f dp m x {ps) f ( v , . 

Ps = 2^9xA x j—^— = f /x(p)- (12) 

x J W \/p 2 + m\{p s ) 



The energy (|TT[) alone gives rise to single-particle energies 



ex(p,Ps) = yp 2 +m x {p s ). (13) 

We take 

u ® = 1 7m- ■ ■ < 14 > 

with £ = p s /po, and a, b, and v adjusted to produce average nuclear ground-state properties. 
The role of the denominator in ( |1~4"D is to prevent supraluminous behavior at high densities. 
To the energies (|13"D, we add in the system frame the gradient, isospin, and Coulomb cor- 
rections, that contribute to the forces, but drop out from collision integrals and velocities, 

e x = lx + A x Ux + t 3X U T + Z x $ , (15) 

where U\ = —a x V 2 (p/po), Ut = a T pr/po, and $ is the Coulomb potential. 

To determine a, b, and v, we required the energy per nucleon to minimize in nuclear 
matter at p = po at the value of e/ p — m^ pa —16.0 MeV for incompressibility K = 210 MeV, 
and at —17.0 MeV for K = 380 MeV. For the higher K, the energetic cost for the surface is 
higher. That leads to difficulties, for the TF theory, in reproducing the average dependence of 
nuclear binding energy on mass number (especially in the low mass region), which we partly 
compensate for with the stronger binding in the infinite- volume limit .Q The parameter sets 
resulting from adjustments are: a = 187.24 MeV, b = 102.62 MeV, and v = 1.6339 for 
K = 210 MeV, and a = 121.26 MeV, b = 52.10 MeV, and v = 2.4624 for K = 380 MeV. 
Generally, the reproducing of the binding-energy curve matters for assessing the excitation 
energy of transients formed in low-energy or peripheral reactions. 



1 A more thorough discussion can be found in [20| |; adjustments of a\ in (0) cannot be done 
without worsening the TF description of measured rms radii. The fact that we employ the average 
of proton and neutron masses for nucleons is of no importance for the issue. 
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In the case of MFs dependent on momentum in their nonrelativistic reduction, we 



also 0)' J = J2x9xA x J j^fx v x = 0, with 



parametrize the energy density in the local frame where baryon flux vanishes [10, 15] (see 

dp 



x 



(16) 



where U is of the form expressed by Eq. fli~4]) , with £ = p/po, and the local particle velocity 
v* x depends on (kinematic) momentum and density through 



P 



I p 2 + m x 



1 I r rnN_ A x g 

1 ^ C m x {l+Xp^/m x y 



(17) 



The energy fll6"|) alone yields the local single-particle energies 



£x(p, p) = m x + / dp'v* x + A x 
Jo 



where 



dp 



(18) 



(19) 



In their nonrelativistic reduction, the energies ([18]) are similar to the energies proposed for 
the nonrelativistic transport by Bertsch, Das Gupta et al. [E21E3. Principally, for matter 



without local reflection symmetry in momentum space, there is a correction term in Eq. (|18|) 
from the condition J = 0, which we ignore as it would seriously complicate our calculation. 
Suprisingly, that term is never mentioned in the context of nonrelativistic calculations. 
On the other hand, in practice, the omission of that term causes no problems with the 
energy-momentum conservation of any relevance, as will be indicated. Different sets of 
parameters for (|16|) , giving different values for the group effective mass |24| in normal matter 



at the Fermi surface (m* = p F /v F ) and for the incompressibility, are exhibited in Table |. 



C. The Thomas-Fermi Equations 



The requirement that the energy (0) is minimal in the ground state, for a nucleus with 
a definite number of protons and neutrons, yields the set of TF equations: 

= ~e P (P F (P P )) ~ a, V 2 (£-) + ^ + * - A* , (20) 

VPo/ 4 po 

= In (p F ( Pn )) - ai V 2 U) - ^ - p n , (21) 

VPo/ 4 po 

and the condition V„p = at the edge of the density distribution. In the equations, p p and 
p n are the Lagrange multipliers for the proton and neutron numbers, respectively. 
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The role of the derivative correction in (0), fl21|), and flT5|) is to reduce the effect of 
the negative MF when the density distribution in the vicinity is primarily concave and to 
enhance the effect of the field when the density is convex. Such a result would be obtained for 
a finite-range effective two-body interaction convoluted with density expanded in position to 
second order. Not surprisingly, the derivative correction is small but it becomes important 
when the energies balance, permitting an adequate description of the density in the ground 
state. For nucleons with high momenta relative to a system, the sign of the correction, 
obviously, becomes an issue. 

In finding the density profile, it is convenient to transform the TF equations into: 

^^^ F + * + i »-^J- < 22 > 

~F i ^ ~F , a T Pp — Pn /OQ x 

Pp - Pn = t + $ - e n + -7T — n • (23) 

1 Po 

The net density profile p(r) may be obtained by starting Eq. ([£3) with some density at r = 0. 
At any r, separate p p and p n may be found from fl23|) and <3> can be obtained from the Gauss' 
law. The acceptable starting p(r = 0) is the one for which p = is reached in the solution 
simultaneously with dp/dr = 0. The chemical potentials are adjusted until the required 
proton and neutron numbers are obtained. At the end, $ and p v may be renormalized, so 
that $ — ► as r —>■ oo. 

Figure [l] shows the calculated proton and neutron density profiles for a moderate and 
large nucleus and MFs corresponding to K = 210 MeV, with and without momentum- 
dependence. Also the empirical charge densities are shown. 



D. Optical Potential 

Comparing relativistic MFs, e.g. when employing both scalar and vector potentials with 
momentum dependence, and nonrelativistic fields, poses some difficulty. Wave approaches 
are summarized at times in terms of the Schrodinger equivalent potential. Feldmeier and 
Lindner [^£| (but see also [p7)l ) suggested the definition of the optical potential in nuclear 



matter as simply the difference: 

U opt (p) = e(p) - m - T{p) , (24) 

where T is kinetic energy at the same momentum p in free space. The momentum can be 
expressed in terms of net energy (such as that of incident nucleon), whether in a relativis- 
tic or nonrelativistic approach, yielding U opt (e). The resulting potential is similar to the 
Schrodinger-equivalent potential at low to moderate momenta and has a more satisfactory 
high-energy limit. 

Feldmeier and Lindner used fl2~4|) to represent the nonrelativistic results in the work by 
Perey and Perey p8 in a combination with the results of an analysis of proton scattering 



data, in terms of relativistic MF potentials, by Hama et al. ||. The best fit |26| with an 
analytic ansatz is shown by a solid line in Fig. [| The dashed lines in the figure show 
potentials in normal matter for three of our K = 210 MeV parametrizations from Table |. 
The omitted K = 380 MeV and m*/m = 0.70 potential is practically indistinguishable from 
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the displayed K = 210 MeV and m*/m = 0.70 potential. This is because the dependence 
of the velocity on momentum and density is forced to be the same for the two potentials, 
cf. Eq. ([T7|) and Table |. In Fig. 0, the linear dependence of the potentials on energy in the 
low-energy region indicates that the effective mass provides a good characterization of those 
potentials in that region. The use of scalar (rather than vector) density in MFs without 
a momentum dependence in the nonrelativistic reduction gives rise to a weak dependence 
(m*/m = 0.98 at e F ) for relativistic energy, as illustrated by the dotted line in Fig. 

With regard to other than normal densities, in some analysis of high-energy scattering, 
see e.g. [Q, it was found that potentials of a bottle shape described the data better than stan- 
dard Woods-Saxon potentials, with the interior at times repulsive while the surface region 
attractive. (That was, in fact, obtained in the analysis || of Hama et al. but, primarily, as 
a byproduct since the parametrization did not have enough flexibility to address the surface 
directly.) Such results indicate differences in the momentum dependence at different densi- 
ties, with the potential changing sign at different energies at the different densities. That 
contradicts a naive expectation that the potential at a lower density may be obtained by 
rescaling of the potential at the normal density; the result is actually predicted by nuclear- 
matter calculations. Figure |3] displays the potentials from early Brueckner-Hartree-Fock [29 



and variational P0[ calculations at p and po/2. It is seen that the p /2 results continue to be 



attractive up to higher energies than the po results. That is also present in our parametriza- 
tions. The increasing slope of a potential with energy, as the density rises, generally indicates 
the increasing velocities with the rise in density, on account of the momentum dependence. 

The flow data, that we shall analyze, will turn out to be sensitive to the potential at 
supranormal densities. After the analysis, we will compare our results to the microscopic 
calculations done at such densities. Feldmeier and Lindner attribute an error of the order 
of 5 MeV to their optical potential fit. In our analysis, we will try not to be guided directly 
by the nucleon data, but rather explore what the flow data alone may tell us about the 
potential. 

E. Lattice Hamiltonian 

Advancing of physical conclusions from comparisons to data depends on advances in the 
methods. In solving the Boltzmann equation (Jl]), the phase-space density is represented in 
terms of a set of (^-functions, or test-particles, 

/(p, r,i) = i^ *(p - p k (t)) 5(r - r k (t)) , (25) 

k 

where M is the number of test-particles per particle. Such a representation is possible when, 
in its use, the phase-space density gets integrated over phase-space volumes incorporating 
many test-particles. In the solution of (|1|), the changes in the density with time due to 
collisions are accounted for by a Monte-Carlo procedure, while the changes described by 
the l.h.s. of (|I|) are accounted for by requiring that the test-particles follow the Hamilton's 
equations: 

1 — (26) 

dt dp k ' 
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d Pk = (27) 
dt dr k ' 1 ' 

For the Hamilton's equations, an averaging in space must be done to obtain single- 
particle energies and gradients must be calculated. Often, a mesh in space is employed. 
The representation of the phase-space distribution in terms of the test-particles leads to 
statistical fluctuations in the driving terms in the Hamilton's equations. Due to these 
fluctuations the net energy and net momentum of the system fluctuate. As a convex function 
of momenta, the energy grows with time as an effect of the diffusive process. The problem 
is more serious for momentum-dependent than momentum-independent MFs because of 
stronger fluctuations. The former fields depend both on position and momentum of the test- 
particles and, additionally, the driving terms fluctuate in two Hamilton's equations rather 
than only in one. The problem may be so severe that, for momentum-dependent fields and 
for ad hoc methods of calculating the single-particle energies and gradients, in a slow central 
reaction at few tens of MeV/nucleon, the spurious energy gain may get close to the net 
excitation energy. Another example when fluctuations can be very harmful, is when using 
the MF to simulate the chiral phase transition. If such transition is not simulated though, 
the situation generally quite improves with the rise in the energy of a simulated reaction. 

The fluctuations may be reduced by increasing the number of test-particles Af, for the 
sake of improving the accuracy of a calculation, but the inverse square-root covergence is 
slow. That is worsened by the fact that, otherwise, an improvement in the accuracy requires 
the reduction in the spatial region for averaging and gradients. Thus, other ways of reducing 
fluctuations and their consequences and of accelerating the convergence of calculations must 
be sought, on top of using a large N '. The lattice Hamiltonian method proposed first by 
Lenk and Pandharipande for momentum-independent fields accomplishes three things. 
The method ensures a continuous variation with position of the test-particle contributions to 
local densities, a continuous variation of the test-particle momenta with time, and it imposes 
a net energy constraint on the dynamics. The constraint plays a local role, on the scale of 
the discretization in space; this is achieved through a consistency between the computation 
of the test-particle contributions to the densities and of the forces on the particles. As 
a rule, when following the method, an excellant momentum conservation also is observed 
for isolated fragments. Here, we generalize the method to the momentum-dependent and 
relativistic fields. The generalized method was already utilized in \TE\ and was in parallel 
developed in plf. 



Within the computational region, we introduce a mesh with nodes at r a separated by 
Ali, i = 1,2,3, in three carthesian directions. With each of the nodes, we associate a form 
factor < S a < 1, continuous and piecewise differentiable within the computational region 
and concentrated around r a . We require the form factors to satisfy: 



J2S a (r-r a ) = l, (28) 



for all r within the region, and 



drS a (r) = AV , (29) 
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where AV = Al% AI2 A/3. The second requirement is that, for averaging, every node gets its 
share of the volume and the first requirement is that every test-particle within the region is 
fully accounted for. The average Wigner function for a node, in terms of S a , is 

7(P,r a )= * ^26(p-p k )S a (T k -T a ). (30) 



k 



The approximate energy (lattice Hamiltonian) in terms of the spatially averaged Wigner 
functions is then 

E = AV (H7(r a )} + ^Pt + ^ pM *(r a ) 



Q 

a 1 



4p 

where $ satisfies a discretized Poisson equation |32] with p. 

The single-particle energy for a test-particle is the variation of E with respect to the 
particle number. The energy turns out to be a weighted average of energies associated with 
the nodes in the vicinity: 

e(P>r) = ^S Q (r - r a )e(p,r Q ) , (32) 

a 

with 

e x (p,r Q ) = i x (p,{J(r a )}) + A x U 1 (r a )+t 3X U T (r a ) + Z x $(r a ), (33) 
where Ut{?o) = cltPt^o)/ Po an d 

Ui (r a ) = — ( 2p(r a ) - p(r a + Ah u { ) - p(r Q - Al t u*)) . (34) 

po ^ {Ali} 2 \ ) 

The derivatives of the the single-particle energy ( |32|) yield an expression for the velocity, as 
an analogous average to that for the energy, 

v(p,r) = ^S Q (r - r Q )v(p,r Q ) , (35) 

a, 

and an expression for the force, amounting to a prescription for the gradient, 

<9e(p,r) ^dS a (r-r a )_ 
~ -flF" = ^ * £(P ' * a) • (36) 

a 

In the simulations, we take S(r) = Y[i9i( r i/^h)- I n the interior of our computational 
area, we use g(x) = 0.5 for \x\ < 0.5, g(x) = 0.75 — 0.5|x| for 0.5 < \x\ < 1.5, and g{x) = 
for 1.5 < \x\. At a forward edge, we use g(x) = 0.75 + 0.5x for —1.5 < x < 0.5, and g(x) = 
outside of that interval. The accuracy of energy conservation for momentum-dependent 
fields, following the lattice Hamiltonian method, is illustrated in Table |T|. 
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F. Collision Rates 



Different expectations might be held with regard to in-medium scatterings. For pertur- 
bative processes, matrix elements in the cross section would be expected to be similar in 
medium to those in free-space. However, for strong interactions, with a number of partial 
waves participating, the cross section might be invariant as a geometric quantity, if not 
for the fact that the space for the independent scatterings might get limited, in a dense 
medium. We explore two types of cross sections in the medium: the free cross sections and 
the cross-sections reduced in such a manner that their radii are limited by the interparticle 
distance, 

a = (i tanh(a /ree / a ) , (37) 

where a = p~ 2 / 3 . The reduction is only applied to the elastic cross sections, in order to 
maintain the detailed balance relations in the medium. Some support for the reduced cross 
sections stems from measurements of the linear momentum transfer and of the ERAT cross 



sections ||35|| . It may be mentioned, regarding the self-consistent nuclear- matter calculations, 
that a reliance on the quasiparticle approximation excludes, in practice, the ability to see 
any changes in the scattering rates due to the overlap of scattering regions; it is necessary 



to incorporate a spreading of the states p3| , |34| . 

The momentum dependence of the MFs can make the calculation of the collision rates 
difficult because of the deformation of the energy shell. We adopt a simplifying approxima- 
tion of de/dT « const, where e and T (kinetic energy) pertain to the same momentum, for 
collisions that locally contribute mostly to the transport. Our approximation is a relativis- 
ts analog of the effective mass approximation. Its validity, throughout the energy range, 
is suggested by the slowing of the variation in the MF and in the optical potential as the 
energy increases, cf. Eq. (^) and Fig. || Within the approximation, we can calculate the 
rates following the same procedure as for momentum- independent fields |!(J (see also [|37j1), 
colliding test-particles within the cells between the nodes, renormalizing only the rates for 
the average change in velocities. One way of assessing the quality of the approximation is 
by examining the accuracy of energy conservation. It is seen in Table that the error does 
not exceed that from integrating the MF motion; in some cases the two errors appear to 
compensate. 



III. ELLIPTIC FLOW 

A. Flow and the Momentum Dependence of the Mean Field 

At few hundred MeV/nucleon, the best probes for the average forces in heavy ion colisions 
are at present the directional features of the transverse collective motion. The time for the 
development of the directional features is limited by the passage of spectators near the 
reaction zone. The overall strength of the motion is not such a good probe because the 
time for the development of the motion as a whole is not comparably limited. Weaker 
forces in the matter are likely to be correlated with longer expansion times and stronger 
forces with shorter times, with the difference in the times compensating for the difference 
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in the forces in the net strength of motion. Besides the size of the acting forces, an issue 
in the reactions is the nature of these forces. The forces may be associated with MF, 
but the matter can also be pushed by thermal pressure. Varying contributions were found 
in early central-reaction simulations |38| utlilizing different combinations of density and 



momentum dependence of the MF. A weak momentum dependence allows for a greater 
increase in density and a better equilibration and, hence, for a greater role for the thermal 
pressure. In assessing the momentum dependence, we shall try to circumvent the question 
of equilibration in central collisions and to follow up on that question in the light of the new 
results, as well as to reexamine EOS, in the future. 

The flow anisotropics may be quantified in terms of the average transverse momentum 
component in the reaction plane or in terms of the eigenvalues of the transverse momentum 
tensor, as functions of rapidity. More flexibility in the quantification give the moments of 
the azimuthal angle relative to the reaction plane, 

v n = (cos (n<f>)) , (38) 

as these can be investigated both as functions of rapidity and transverse momentum. 
At midrapidity, away from spectator contributions, the lowest nonvanishing moment in 
symmetric systems is v 2 describing ellipticity of the particle distribution; the derivative of 
t>i, dvi/dy, may be studied in that region, too. 

Figure f| shows the dependence of the midrapidity v% on impact parameter, obtained in 
the Boltzmann-equation simulations of Au + Au collisions at 400 MeV/nucleon. Mean fields 
with and without momentum dependence, corresponding to different incompressibilities, 
were utilized in those simulations. The negative values of Vi for all simulations indicate 
a preference for particle emission out of the reaction plane, towards 90° and 270°. For low 
impact parameters, the system approaches azimuthal symmetry in space and the magnitude 
of v 2 decreases. In semicentral collisions, b < 6 fm, the values of v-i from simulations with 
different MFs cluster differently at different impact parameters, reflecting an interplay of the 
geometry, of the dependence of MF on density and momentum, and of the thermal pressure. 
Beyond that lower-impact parameter region, however, the results of simulations utilizing MF 
with and without momentum dependence, exhibiting distinctly dissimilar behavior, clearly 
separate. Without momentum dependence in the MF, the values of v% first stabilize and then 
drop in magnitude. With the momentum dependence, the values of V2 continue to increase 
in magnitude up to very high impact parameters, in a roughly linear fashion. These results 
give rise to an expectation that the measurements of ellipticity at high impact parameters 
can be used to test the momentum dependence of the MF in collisions. 

Due to its very nature, one might further hope that the momentum dependence of 
MF would be revealed in the features of emission of particles with high momenta. In the 
past, in fact, we have indicated fK| that the behavior of the first-order flow with transverse 
momentum, Vi(p ), could be used to discern the momentum-dependent from the momentum- 
independent MFs in collisions. To our knowledge, though, no experimental analysis of those 
coefficients has been done within the energy range where the degrees of freedom, playing 
a role in the collisions, are still under some control. However, the intermediate-energy data 
on the dependence of ellipticity on p x exist p9|,^0| and, following the discussion, one might 



hope to access the momentum dependence of the MF in v 2 at high p and b. 
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B. Data Comparisons 



Figure |5] compares the p -1 - dependence of the measured and calculated ratios of out-of- 
plane to in-plane midrapidity proton-yields (0.35 < y/ybeam < 0.65), Rn = (iV(90°) + 
iV(-90°))/(iV(0 o ) + iV(180°)), in 400 MeV/nucleon 209 Bi + 209 Bi collisions at b ~ 8.7 fm, 
as computed from the lowest Fourier coefficients, Rn = (1 — i>2)/(l + 1*2)- We utilize the 
KaoS [4(J data in our comparison, rather than the nominally similar LAND [39] data, 



because of the wider available impact-parameter range and because of an absolute normal- 
ization of the KaoS data. The top panel of Fig. || compares the data (squares) to the cal- 
culations (lines) done for in-medium cross sections and different MFs which all correspond 
to the incompressibility of K = 210 MeV. The experimental ratio rises rapidly with the 
transverse momentum above 300 MeV/c, reaching values higher than 2 above 500 MeV/c. 
The ratio from simulations without momentum dependence in the MF stays relatively flat 
with values below 1.4 up to momenta of 700 MeV/c. It may be mentioned that the continu- 
ity of the momentum distribution enforces |41| a quadratic behavior of v 2 , and thus of Rn, 
around p = 0. With the momentum dependence in the MF, a strong sensitivity to the 
details in that dependence is observed in the p^-dependence of Rn, with the KaoS data at 
the displayed momenta favoring the MF parametrization characterized by m*/m = 0.70. 
The validity of conclusions on the momentum dependence hinges on the sensitivity of the 
results to the incompressibility and to cross sections. The sensitivities are tested in the bot- 
tom panel of Fig. []. Simulations with MFs without momentum dependence, corresponding 
to different incompressibilities, and the cascade model yield all practically the same results 
for RnIp' 1 ) at high p- 1 . A sensitivity to the in-medium NN cross sections is observed, but it 
is weak in the high-momentum region. 

We next examine whether similar conclusions can be drawn at other beam energies. Thus, 
the top panel of Fig. ^] compares the p^-dependence of the measured and calculated midra- 
pidity ratios in 700 MeV/nucleon 209 Bi + 209 Bi collisions at b ~ 8.6 fm. Again, the ratios 
from simulations with momentum-independent MFs are quite incompatible with the data. 
As to the details in the momentum dependence, the region of high sensitivity to those de- 
tails shifts towards higher transverse momenta, compared to collisions at 400 MeV/nucleon. 
We should indicate that our ability to address very high transverse momenta is limited by 
the finite statistics in the simulations. The data, nonetheless, eliminate such weak mo- 
mentum dependence as that represented by the m*/m = 0.79 parametrization. Further- 
more, the 700 MeV/nucleon data marginally favor the parametrization characterized by 
m*/m = 0.65 over m*/m = 0.70. It is seen in Fig. |6|, that some sensitivity to just the 
density dependence of the MF develops at 700 MeV/nucleon in the midperipheral collisions. 

Some further shift of the region in p 1 - where ellipticity strongly discriminates different 
momentum-dependences of the MF is observed in the 1000 MeV/nucleon Bi + Bi collisions, 
see the top panel of Fig. |7|. While the predictions for the momentum-dependent MFs clearly 
separate from those for the momentum-independent MFs and are, by far, more consistent 
with the data HOI, the separation between the predictions for the m*/m < 0.70 and m*/m = 
0.79 MFs is not that large at p 1 < 1000 MeV/c. 

Each time, so far, we compared the theory to data at about the same impact parameter 
at the three beam energies. The bottom panel of Fig. [7| illustrates the quality of data 
description, by the calculation with m*/m = 0.70, at different impact parameters. We should 
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mention that we refrain from deciding on the momentum dependence of the MF based on 



comparisons to the highest available impact parameters |^0[ , as nonlinearities in the variation 
of quantities under the averaging over impact parameter in the marginal region may affect 
the results. We should also mention that for some combinations of impact parameter and 



beam energy, the values of out-of- to in-plane ratios Rn from the KaoS measurements jllj fall 
below 1 at low jr 1 . We do not observe such a behavior in the calculations. The preliminary 
FOPI data on the p^-dependence of elliptic flow do not exhibit such a behavior, either p2j . 



Summarizing the analysis of simulations and the comparisons to data thus far, we find 
that the ellipticity at high impact parameters and transverse momenta is a very sensitive 
probe of the momentum dependence of the nucleonic MF. The large measured values of 
out-of-plane to in-plane ratios (or ellipticity) at midrapidity cannot be explained without 
the momentum-dependence. Selection of data favors an MF parametrization characterized 
by m*/m ~ 0.70. Coincidentally, such parametrization agrees with the information on 
the nucleonic MF at po from nucleon scattering. A natural question to ask at this stage is 
whether the measurements of the ellipticity represent just another way of accessing the same 
information or whether these measurements have the potential to expand on the information 
from nucleon scattering. 

A surmise regarding the same information as from scattering might be based on the fact 
that the densities reached in peripheral collisions are lower than in the central ||41||. In sim- 



ulations with the momentum-dependent MF at b ~ 8.7 fm, we find, though, the maximal 
densities of p max /po ~ 1-85, 2.20, and 2.40, at 400, 700, and 1000 MeV/nucleon, respec- 
tively. However, the spatial volume with such densities might not be large or time span for 
the densities long enough to affect the particle emission, or otherwise the maximal densities 
might not be relevant. To assess whether the elliptic flow tests the momentum dependence 
more at supranormal or at normal and subnormal densities, we have carried out the simu- 
lations with an m*/m = 0.70 and K = 210 MeV MF modified by making the momentum 
dependence at p > p follow the dependence at p . The modification was accomplished by 
demanding that the velocity in the matter in ( |i~6|) ceased to change (i.e. increase) beyond 
the normal density: 

v* X (P,0=v*x(p,1), for£ = p/po>l- (39) 

As the freezing of the momentum dependence softens the zero-temperature energy density 
above p , we compensated that softening by adding a repulsive term to the potential U 
in (|16D and ([14]) , of the form 25 MeV (^ — l) a35 , for £ > 1, roughly restoring the density 
dependence of the energy from before the modification. Note that freezing the momentum 
dependence at zero rather than normal density would yield a momentum-independent MF. 
The results of the simulations with the momentum dependence of MF frozen above po are 
shown in the bottom panel of Fig. |6|. It is seen that the yield ratios are now well 
below data at high transverse momenta. The results are, in fact, a bit closer to the results 
obtained with the momentum-independent MF than to the results obtained with the MF 
with the standard dependence. Clearly, the elliptic flow tests the momentum dependence at 
supranormal densities. 

In the view of the last finding, the preceeding comparisons to data at 700 and 
1000 MeV/nucleon and at high transverse momenta suggest that the momentum depen- 
dence in the so-far favored m*/m = 0.70 MF parametrization may not actually grow fast 
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enough with the density at the highest reached densities and momenta. We next investigate 
reasons for the sensitivity to MF at high densities. 

C. Dynamics and Anisotropy at High Transverse Momenta 

To understand the association of high-momentum anisotropics with the momentum de- 
pendence of MF at high densities, we examine the production rate for high-momentum 
midrapidity protons in midperipheral collisions. In the upper panel of Fig. ^ we plot the 
emission rate as a function of time for protons with p 1 - > 550 MeV/c in the K = 210 MeV 
m* /m = 0.70 Bi + Bi simulations at 400 MeV/nucleon and b = 8.7 fm; we also plot in that 
figure the density at the system center. The emission time is defined as the instant of last 
collision. It is seen that the peak emission is strongly correlated with the system reaching 
a high density at the center. For comparison, we plot in the bottom panel of Fig. |8| the 
emission rate for all midrapidity protons. The emission for all protons only sets in during 
the high-density stage. Thus, at t = 30 fm/c, when 90% of all high-momentum protons 
have already been emitted and central density has fallen to low values, fewer than 30% of all 
midrapidity protons have been emitted. Of all midrapidity protons, these are the high-p- 1 
protons that are most directly emitted from the high-density zone during the overlap of the 
nuclei. For a view of the overall progress of the Bi + Bi reaction, see Fig. [|. 

A possible cause for the large high-p -1 anisotropies obtained with the momentum- 
dependent MFs, larger than with the momentum-independent MFs, might be a stronger cor- 
relation of the high-p^ emission with the overlap of the nuclei for the momentum-dependent 
fields. The spectator pieces are, generally, expected to shadow the emission within the 
reaction plane. To clarify the issue, in addition to the high-p -1 rate for the momentum- 
dependent MF, we plot in the top panel of Fig. ^ the high-p -1 emission rate obtained with 
the K = 380 MeV momentum-independent MF. While there are some statistically signifi- 
cant differences between the rates, their overall time dependencies are fairly similar; in no 
way the small differences could explain the huge differences in anisotropies at high p ± found 
in Fig. |5j. The time dependencies for the central densities and for the density distributions 
are quite similar for the two MFs, difficult to distinguish by eye. Thus, some differences in 
the directionality of the emission for the two MFs must exist at the same emission times, 
and be responsible for the large ratios in Figs. |5|-|7|. 

To more understand the emission, we shall examine v<i = (cos 20), for the high- 
momentum protons, as a function of time. Before that, we need to discuss expectations 
regarding t> 2 i n simple situations. For emission strictly towards 90° and 270° relative to 
the reaction plane, we expect t>2 = —1. For a model emission from rough flat surfaces 
with normal directions towards 90° and 270°, as suggested by the 10-20 fm/c top panels 
m Fig. g, the emission pattern is dN/dcj) oc |sin0|. The value for the coefficient is then 
t> 2 = —1/3. (The roughness scale, that should be much smaller than the surface size, is the 
mean-free-path in our case.) Another interesting value of v<i is that for the most anisotropic 
distribution describable in terms of two Fourier cofficients, i.e. dN/dcf) oc (1 — cos 20), for 
which t>2 = — l/2. 

Figure |Tl] shows the time-dependence of — t>2 for high-p -1 protons from the 
400 MeV/nucleon Bi + Bi simulations at b = 8.7 fm employing the two recently dis- 
cussed MFs. Both the average values (top panel) for protons emitted in the vicinity of 
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a given time, — (cos 20) (t), as well as the values (bottom panel) for all protons emitted up 
to a given time, — (cos 20) t / <4 , are shown. It is apparent from Fig. [| and from the bottom 
panel of Fig. [10] that the values of t>2 for particles emitted past ~ 25 fm/c have already little 
effect on the overall value for all particles. 

As is seen in the top panel of Fig. [10], the emission for the momentum-independent 
MF starts out with values of — t>2 ~ 0.15, in the vicinity of a half of what is expected 
for a flat surface parallel to the reaction plane. It appears that the emission geometry is 
more complicated than in the latter case, with a finite size of roughness for the surface 
and with the surface oriented towards different angles as may be apparent from the bottom 
5-15 fm/c panels in Fig. |9[. In contrast to the momentum-independent MF, the emission for 
the momentum-dependent MF starts out with much higher values of — t>2 ~ 0.50, well above 
the expectation for a flat surface. As geometry is virtually the same as for the first MF, 
it must be a difference in the field within which the particles move that plays a focussing 
role. Basically, the effect may be understood in terms of a relatively strong repulsive MF 
felt by the high-momentum particles in the second case in the high-density zone. As the 
particles move out, they feel a strong gradient towards the normal of the emitting surface. 
This gradient pushes out the otherwise rapidly falling distribution of emitted particles in the 
transverse direction. The highest transverse momenta are acquired by particles best aligned 
with the gradient which gives rise to the enhanced anisotropy at the high momenta. As 
the spectator pieces pass by the center of the system, there is a burst of particles that were 
moving in the transverse direction parallel to the reaction in the participant zone and were 
trapped till then. That lowers the instantaneous — t>2 values in Fig. |10] around t ~ 20 fm/c 
for either of the employed MFs. Later, the instantaneous values recover somewhat, possibly 
due to left-over particles from the early stage that just underwent some secondary collisions 
and got reclassified with regard to their emission time. 

Optical potential fields at a high-density stage, in simulations of the Bi + Bi reaction 
employing different MFs are shown in Fig. [II] for high-p^ nucleons. The MFs include the 
momentum-independent MF characterized by K = 380 MeV, the MF characterized by 
m*/m = 0.70 and K = 210 MeV MF, and that last MF with the momentum dependence 
frozen above po- The MF with frozen momentum dependence fails to describe the magni- 
tude of elliptic anisotropy at high p in the semiperipheral Bi + Bi at 700 MeV/nucleon 
and, likewise, at 400 MeV/nucleon. In the high-density zone, the optical potential for 
p 1 - = 600 MeV/c nucleons is strongly repulsive for the standard MF, reaching values of 
~ 55 MeV. Still higher values are reached for higher momenta. Notably, in normal mat- 
ter, the optical potential for that MF, as well as the potential from nucleon scattering, 
vanishes in the vicinity of 600 MeV/c (see also Table |). In the case of the frozen momen- 
tum dependence, the optical potential in the high-density zone reaches only ~ 15 MeV. 
Finally, for the momentum-independent MF, the optical potential stays mildly attractive 
at 400 MeV/nucleon in the high-density zone, reaching maximal values there in the vicin- 
ity of ~ —5 MeV. The attractive potential could principally play a defocussing role in the 
emission; the net results for anisotropy are, nonetheless, similar to the results without MF, 
as was demonstrated earlier. The buldging out of the momentum distribution, due to a re- 
pulsive field for high momenta, suggests the possibility for positive values of the fourth 
Fourier-coefficient at midrapidity, V4 = (cos 40); the coefficient would account for a fine 
structure in the enhancement around 90° and 270°. The coefficient is shown as a function 
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of transverse momentum in Fig. [12] for the m* jm = 0.70 K = 210 MeV MF and for the 
momentum-independent K = 380 MeV MF. For the momentum-dependent MF the coeffi- 
cient rises indeed rather rapidly with momentum, whereas it remains close to zero for the 
momentum-independent MF. 

One outstanding question, which we need to address, is whether we can actually draw 
conclusions on the effective mass in ground-state matter, be that in an average sense. The ef- 
fective mass was, so far, quite useful in labeling the different MFs. To address the question, 
we construct an MF parametrization requiring that the effective mass at Fermi energy in 
the ground-state matter is m* jm = 0.65 and that the optical potential vanishes at the same 
momentum as for our previous to*/ m = 0.70 MF, i.e. at Pq ~ 680 MeV/ c, cf. Table |. The re- 
quirement makes the potentials appear similar at intermediate momenta; at high momenta, 
though, the m* jm = 0.65 potential becomes softer than m* jm = 0.70, approaching a lower 
asymptotic value, see the table. We now turn to the predictions for ellipticity. Figure |13 
shows the calculated ratio of in to out-of plane proton emission at midrapidity as a function 
of transverse momentum in midperipheral Bi + Bi reactions at 400 and 700 MeV/nucleon, 
for the new and previous MFs, together with the data. At 400 MeV/nucleon, the predic- 
tions obtained using the m* jm = 0.70 MF and the new m*jm = 0.65 MF are fairly similar. 
With the increase in beam energy, however, the predictions separate, with these for the new 
m* jm = 0.65 MF falling below those for the m* jm = 0.70 MF and below the data. Clearly, 
at the considered energies, the high-momentum elliptic flow is sensitive to the momentum- 
dependence of MF at intermediate and high momenta but not at momenta as low as the 
ground-state Fermi momentum. On the other hand, given the requirements of the energy 
and density for the ground state, and of the optical potential becoming repulsive but not 
too repulsive with the increase in momentum to yield right anisotropies, it is difficult to 
put forward a MF with m* outside of the interval m* jm ~ (0.63 — 0.73), unless that MF 
were to have rather discontinuous behavior. Of course, on account of the dispersion rela- 
tion involving single-particle energies and scattering rates and the discontinuity in the Pauli 
blocking at zero temperature, there might be variations in m* in the very vicinity of the 
Fermi momentum, but these would be quite local. Now we turn to a direct comparison of 
our results to microscopic calculations. 



IV. COMPARISONS TO MICROSCOPIC OPTICAL-POTENTIAL 

CALCULATIONS 

The best way to test the momentum-dependence in microscopic MFs would be to imple- 
ment these MFs into a transport calculation and confront the elliptic-flow results with data. 
Given though the preliminary nature of our investigation, where the elliptic flow is exploited 
for the first time in exploring the momentum dependence, we will merely compare directly 
the utilized optical potentials to those calculated in the literature, at densities reached in 
reactions and at momenta where anisotropies are observed. That will also give the opportu- 
nity to take a closer look at the potentials in simulations, presented so far only for normal 
and subnormal densities. We shall comment on some difficulties in the comparisons on the 
way. 
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A. Optical Potentials from Dirac-Brueckner-Hartree-Fock Calculations 



The first microscopic optical potentials that we consider are those obtained within the 
Dirac-Brueckner-Hartree-Fock (DBHF) approach [p7|,|4"3|. In that approach, an in- medium 



Thompson equation is solved with a realistic relativistic interaction. The equation accounts 
for the effects of Pauli principle and of self-consistent scalar and vector MFs. Rise in the 
scalar MF with density leads to a weakening of the scalar exchange in the interaction and 
contributes to the nuclear saturation E3|. In Fig. [TJ], we compare the DBHF optical po- 



tentials obtained using the Bonn-A interaction to our parametrizations. The potentials 



from a relativistic approach may be compared as functions of momentum. The short-dashed 
lines in the figure represent the potential determined in [0,fE| for a broad range of densities 



assuming momentum-independent scalar and vector MFs. The short-dashed-dotted lines 
represent the optical potential determined in || for a more narrow range of densities as- 
suming a parametrized form of the momentum dependence for the scalar and vector MFs. 
The momentum-dependence of the optical potential for the momentum-independent MFs 
stems from a large magnitude of the two MFs that largely cancel out while differently en- 
tering the single-particle energies. That effect also largely contributes to the momentum 
dependence of the optical potential in the case of momentum-dependent MFs; see the simi- 
larity of the two potentials in Fig. [TJ]. The solid and long-dashed lines in the figure represent, 



respectively, our m*/m = 0.70 and first m*/m = 0.65 K = 210 MeV parametrizations. 



It may be seen in Fig. ITJI that, at densities and high-momenta explored in the 
400 MeV/nucleon collisions, the DBHF potential || is fairly close to that from our 
m* /m = 0.70 parametrization. The DBHF potential f7| would likely give a too large high-p- 1 
elliptic anisotropy at 400 MeV/nucleon. As the density increases, the DBHF potential |7| 
rises more and more above our m*/m = 0.70 potential. However, at the higher two beam 
energies in Bi + Bi collisions, involving somewhat higher densities, we did not explore how 
strong the MF would need to be to yield excessive anisotropies. Thus, we cannot classify 
the behavior of the DBHF potential at the increasing densities as unrealistic. We should 
note that the optical potential in our m*/m = 0.70 parametrization appears more repulsive 
at p = 650 MeV/c and at p/po ~ 1.5 in the simulation in Fig. O than at such momen- 
tum at p/po = 1-5 in Fig. [14]. Basically, the effect may be understood in the following way.0 
As particles in a reaction are locally at higher typical momenta than in the zero-temperature 
matter, they feel, on the average, a more repulsive potential. By self-consistency this makes 
the potential more repulsive at any fixed momentum. We expect such dynamic shifts for 
all the potentials. In addition, a given fixed momentum in the overall frame corresponds to 
different momenta in local frames in a reaction which is similar for different potentials in 
simulations. 

One may notice in Fig. [14] that, at p = po = 0.16 fm -3 and momenta p < 400 MeV/c, 
the DBHF potentials are always lower than the potentials from our parametrizations. This 
might be due to a problem in the DBHF calculations with the thermodynamic consistency 
relation (Hugenholtz-van Hove theorem): 



J This type of an effect was discussed before, in different terms, by Wolter et at. [44 
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in the ground state. Our optical potentials cross at the Fermi momentum in Fig. |14 
required by (f40|). The fact is, however, that the DBHF potentials are a bit too attractive 
at low momenta and the nuclear matter within the DBHF approach saturates at a slightly 
excessive density of 0.18 fm~ 3 p3|P|. Overall, though, the proximity of the DBHF potentials, 
especially for the momentum-dependent MFs, to the parametrization which adequately de- 
scribes the elliptic-flow anisotropies, is rather startling in the tested range of densities and 
momenta. 



B. Potential by Baldo et al. 

Next potential, being considered, is that obtained within the nonrelativistic Brueckner- 
Bethe-Goldstone (BBG) approach with the Paris interaction In that approach, to the 
Brueckner-Hartree-Fock term in the optical potential a so-called rearrangement term is 
added which accounts for the influence of the particle, for which the potential is evaluated, 
on the correlations between other particles in the medium. At low momenta, the rear- 
rangement term weakens the dependence of the optical potential on momentum in quite 
an essential manner. The BBG potential, represented by a short-dashed line, is compared 



to our parametrizations as a function of energy at several densities in Fig. [15]. It is seen that, 
at all the supranormal densities in the figure and, essentially, at all energies, the BBG po- 
tential is below our second m* jm = 0.65 parametrization and even below the m*/m = 0.70 
parametrization with the momentum dependence frozen at the higher densities; the poten- 
tials from the two parametrizations are represented by the dotted and long-dash-double- 
dotted lines, respectively. Notably, the nuclear matter in the BBG approach saturates at 
a very high density of p = 0.24 fm~ 3 ||. From our investigations, we can conclude that the 
BBG potential is unacceptably weak at the densities p > 1.5po an d higher energies and mo- 
menta; that potential would yield far too small elliptic anisotropies in the Bi + Bi collisions 
at 400 and 700 MeV/nucleon. 

C. Potentials from Variational Calculations 

The final set of optical potentials, which we examine, is that obtained within the varia- 
tional method for nuclear matter for three combinations of two- and three-nucleon interac- 
tions ||. These combinations included the Friedman and Pandharipande |3(| combination 



of the Urbana t>i4 (UV14) with a three-nucleon interaction (TNI), and the combinations 
of either the Urbana Vu or the Argonne Vu (AV14) nucleon-nucleon interaction with the 
Urbana model VII (UVII) three-nucleon interaction. As criticized in |Bj , the rearrangement 
contribution was not included in determining the optical potentials in ||; such a criticism 
might be also raised against the DBHF calculations. The nuclear matter was constrained 
in the calculations || to saturate at p = 0.157 fm -3 . Earlier Friedman and Pandharipande 
results (UV14+TNI) within the approach, for subnormal densities, were already displayed 
in Fig. |. 
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Figure [16] compares, as a function of energy and at several densities, the optical po- 
tential U for the AV14 + UVII interactions to the potentials for our K = 210 MeV MF 
parametrizations. It is seen that, at the densities in the figure and at high momenta, 
the AV14 + UVII potential, represented by a short-dashed line, is always below the po- 
tential from our second m*/m = 0.65 parametrization; at p ~ 1.7po, it is even below the 
m*/m = 0.70 parametrization with the momentum-dependence frozen. This is then similar 
to the case of the BBG potential. The AV14 + UVII potential is unacceptably weak at 
supranormal densities. Use of this potential in simulations would lead to a serious underes- 
timation of high-momentum elliptic anisotropics in the semiperipheral Bi + Bi collisions at 
400 and at 700 MeV/nucleon. 

Finally, we turn to UV14 potentials. Both the UV14 + TNI and UV14 + UVII po- 
tentials || are compared to the potentials from our parametrizations in Fig. [T7|. It is seen 
that at supranormal densities in the figure, the UV14 + TNI potential, represented by 
a short-dashed line, is very close to the potential from our m*/m = 0.70 parametrization, 
represented by a solid line. The UV14 + UVII potential, represented by a solid-dash-dotted 
line, is below the potential from our second m*/m = 0.65 parametrization, represented 
by a dotted line, at p < 2p . Thus, a use of the UV14 + UVII in a simulation should 
lead to some underestimation of the high-p -1 elliptic anisotropy in the Bi + Bi collision at 
700 MeV/nucleon. 



V. DISCUSSION 

To improve the reliability of transport reaction simulations with momentum-dependent 
mean fields, and thus to enhance the validity of conclusions from comparisons to data, 
we have generalized the lattice Hamiltonian method 117] to the case where the reaction 



transport is formulated within the relativistic Landau theory. We have shown (Fig. f|) that 
the elliptic flow at midrapidity exhibits a particularly strong sensitivity to the mean-field 
momentum dependence in midperipheral to peripheral collisions. A relatively weak sensitiv- 
ity was found in these collisions to the incompressibility of nuclear matter. An additionally 
enhanced sensitivity to the mean-field momentum dependence is exhibited by the elliptic 
flow of particles with high transverse-momentum. Variations in the mean field associated 
with the change in the effective mass from m* /m = 0.79 to m*/m = 0.65 change e.g. the 
high-p^ anisotropy from ~ 2 to ~ 3.5 (Fig. |]). Given that the KaoS Collaboration has mea- 
sured [4~0] the p- 1 - dependence of elliptic anisotropics at three beam energies and at different 



impact parameters in Bi + Bi collisions, we went on to assess the mean-field momentum 
dependence from their data and to assess the origin of the large high-p -1 sensitivity to the 
dependence. 

We have shown in Fig. |8| that the high-density stage in midperipheral collisions of heavy 
nuclei is accompanied by a burst of particles with high transverse momenta. These particles 
probe the high density matter in a quite direct and perturbative manner. Thus, no high- 
p 1 - particles are in practice emitted at any other time in a reaction. Remaining particles 
primarily continue along the beam axis and leave the system at different times. Looking at 
statistics in Fig. [8], it is seen that, even amongst the emitted midrapidity particles, those 
with a high p represent a small fraction of all. Thus, the high-p^ particles do not change the 
high-density environment as they leave. The elliptic anisotropy of the high-p -1 particles is 
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sensitive to the sign and to the magnitude of the optical potential felt by these particles and, 
thus, to the momentum dependence of the mean field. For a repulsive potential, the particles 
are speeded up as they roll off the potential step in a transverse direction, escaping into the 
vacuum. 

The large anisotropies observed by the KaoS Collaboration at high momenta in the 
midperipheral collisions at 400-1000 MeV/nucleon (Figs. 0-0) can only be explained in the 
transport simulations when assuming optical potentials with a momentum dependence that 
strengthens as the nuclear density increases beyond the normal density. For momenta where 
the large anisotropies are observed, the potential must be strongly repulsive; the changes in 
the potential as low as 10 MeV yield observable changes in the anisotropy. The particular 
data constraint the potential in simulations at densities from within the region of ~ (1.4 — 
2.3) po and at momenta which correspond to free-space kinetic energies of ~ (100—400) MeV. 

We have compared optical potentials from the MF parametrizations in the transport sim- 
ulations to the potentials from microscopic nuclear-matter calculations. Potentials in the 
simulations allowed for a different quality of description of the measured anisotropies de- 
pending on beam energy and thus densities reached. We found, in the potential comparisons, 
that two of the microscopic potentials, BBG |6|,§] and AV14 + UVII 0, were unacceptably 
weak at supranormal densities and probed energies. On the other hand, the potentials 
from DBHF calculations and the UV14 + TNI potential from the variational calcu- 

lations || turned out to be surprisingly close, in their momentum and density dependence 
within the accessed region of dependence in reactions, to the parametrized potentials giving 
an acceptable description of the data. 

We hope that present and possible future results from flow in peripheral collisions may 
play a similar role to the results from nucleon-nucleus scattering in constrainting microscopic 
theory, but in this case at supranormal densities. A combination of microscopic theories with 
two sets of constraints can make extrapolations to the regions of uncertainty, such as low 
momenta in matter at high density, more trustworthy. Incidentally, in Figs. [T3] and [T7| there 
is quite a degree of convergence in the low-momentum region between the potentials which 
agree in the high-momentum region. 

With regard to possible future investigations of flow, we encountered here ambiguities 
at low transverse momenta of emitted particles that might become clarified. The study 
should be extended to high transverse momenta, which requires increasing statistics within 
simulations. Moreover, investigations could be extended down [|5] and up |4(| in the beam 
energy. The first-order flow, when analyzed in a similar manner to the second-order, is likely 
to contain a comparable amount of information [[T0|JTT] , |4T|JT^ , |48[1 . 

Better constraints on the momentum dependence of the optical potentials should make 
the determinations of the nuclear equation of state more reliable |T0|jrT| . These searches 
should rather be carried at the more central impact parameters (Fig. 0) where the matter 
is better equilibrated and reaches higher densities. One recent example, illustrating the 
importance of clarifying the momentum dependence in that context, is Ref. |4£| where, in 
the simulations for AGS energies, the baryon optical potential is put to zero above a cut- 
off momentum. This should lead to a discontinuity in the potential growing with density 
and likely generate effects similar to a phase transition that one looks for in that energy 
regime [f46| . Notably, at high incident energies, even in midperipheral collisions, the matter 



is dominated by resonances, so investigations of proton flow are likely to mix information 
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on the potential felt by nucleons with that felt by baryon resonances. 
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TABLES 



TABLE I. Parameter values for the momentum-dependent mean fields. The momentum pg is 
the one for which the optical potential vanishes at the normal density, U op t(po, po) = 0, while U^ t 
is the asymptotic value of the potential at po as p — ► oo. The annotated sets have been used in [15]. 



a 

(MeV) 


b 

(MeV) 


;/ 


c 


A 


m* /m 


Po 
(MeV) 


TTOO 

u opt 

(MeV) 


K 

(MeV) 


Ref. 


185.47 


36.291 


1.5391 


0.83889 


1.0890 


0.65 


585 


55 


210 


1 


m 


185.56 


32.139 


1.5706 


0.96131 


2.1376 


0.65 


680 


23 


210 






209.79 


69.757 


1.4623 


0.64570 


0.95460 


0.70 


680 


40 


210 






214.10 


95.004 


1.4733 


0.37948 


0.55394 


0.79 


900 


25 


210 






123.62 


14.653 


2.8906 


0.83578 


1.0739 


0.65 


580 


56 


380 




151 


128.22 


22.602 


2.5873 


0.64570 


0.95460 


0.70 


685 


39 


380 







TABLE II. Accuracy of energy conservation and supplementary information for simulations 
of b = Bi + Bi collisions at different beam energies T^ eam /A, using the momentum-dependent 
MF given by the first of the value sets listed in Table ffl. The first number for A max (E / A) is the 
maximal deviation from the initial cm energy for a simulation with the MF only; the second number 
is for a simulation with the MF and collisions. The simulations utilized All = AI2 = 0.92 fm and 
J\f = 170. The accuracy improves for lighter nuclei and momentum-independent fields (5 times as 
a rule, in the latter case). 



T b eam/A [GeV/nucleon] 0.040 0.400 10.7 

AZ 3 [fm] 092 092 0.40 

At [fm/c] 0.50 0.30 0.080 

t max [fm/c] 300 110 40 

A max (E/A) [MeV/nucleon] 0.9/0.6 3.9/3.5 29/55 
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FIG. 1. Nucleoli density profiles from solving the TF equations for MFs corresponding to 
K = 210 MeV, together with the empirical charge density profiles for 40 Ca and 208 Pb. The solid 
lines represent the empirical profiles from Ref. p5|| . The long- and short-dashed lines represent 
the proton and neutron profiles, respectively, for the momentum-independent field. The long- and 
short-dash-dotted lines represent the proton and neutron profiles, respectively, for the momen- 
tum-dependent field in Table | that yields m* = 0.70 m. 
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FIG. 2. Real part of the optical potential in normal nuclear matter, as a function of nucleon 
energy. The solid line represents a fit p6| to the information from structure and nucleon scatter- 



ing [ 28 , || . The dashed lines represent the potential for different MF parametrizations in Table | 
for K = 210 MeV. The numbers in the figure indicate the ratio m*/m. For m* jm = 0.65 only the 
first set from the table is represented. The dotted line represents the potential for the MF without 
momentum dependence in its nonrelativistic reduction, corresponding to K = 210 MeV. 
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FIG. 3. Optical potential in nuclear matter, as a function of nucleon energy, obtained in 
Brueckner-Hartee-Fock p9i (top) and variational calculations and in our K = 210 MeV and 
m*/m = 0.70 parametrization. Solid and dashed lines show the potential at po and po/2, respec- 
tively. 



29 



CM 
> 



0.20 



0.15 



0.10 



0.00 



t — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — r 



400 MeV/nucleon Au + Au 



• HM 

o $M 
■ H 

° 



8 n 



□ 



♦ $ 

o 

□ 
□ 



□ 



0.00 
0.0 



U=LJ I I I I I I I I I I I I I I I I I I I I I I I L 



m 

I 



2.5 



5.0 



7.5 



10.0 



12.5 



b [fm] 



FIG. 4. Negative of the ellipticity coefficient at midrapidity, as a function of the impact 
parameter, from simulations of Au + Au collisions at 400 MeV/nucleon. The squares and circles 
represent, respectively, the results for MF without and with momentum dependence (m*/m = 0.65, 
of Ref. [15]). The filled symbols are for the incompressibility K = 380 MeV and the open symbols 
are for K = 210 MeV. 
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FIG. 5. Measured [Q (filled squares) and calculated (lines) ratios of out-of-plane to in-plane 
proton yields at midrapidity (0.35 < y/y beam < 0.65), R N = (N(90°) + iV(-90°))/(iV(0 o ) + 
./V(180 o )), as a function of transverse momentum, in 400 MeV/nucleon 209 Bi + 209 Bi collisions 
at b = 8.7 fm. Both the measured and the calculated ratios are computed from the lowest Fourier 
coefficients of the proton distribution, Rn = (1 — V2)/(l + ^2)- The lines in the top panel represent 
results of the simulations done using in-medium cross sections and those MFs for which the optical 
potentials are shown in Fig. ||. The numbers in the panels indicate the magnitude of effective mass 
for the MFs. The bottom panel shows the sensitivity of the results to the variation of cross sections 
and of incompressibility. The long-dashed, solid, and dotted lines repeat respective results from 
the top panel obtained with in-medium cross sections and K = 210 MeV. The long-dash-dotted 
and short-dash-dotted lines represent additional results obtained, respectively, using the momen- 
tum-independent MF corresponding to K = 380 MeV and using no MF at all. The short-dashed 
and long-dash-double-dotted lines represent the additional results obtained for free cross-sections 
and MFs corresponding to m* /m = 0.65 and to m*/m = 0.70, respectively. 
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FIG. 6. Measured [Q (filled squares) and calculated (lines) ratios of out-of-plane to in-plane 
proton yields at midrapidity (0.35 < y/ybeam < 0.65), as a function of transverse momen- 
tum, in 700 MeV/nucleon 209 Bi + 209 Bi collisions at b = 8.6 fm. Both the measured and the 
calculated ratios are computed from the lowest Fourier coefficients of the proton distribution, 
Rn = (1 — f2)/(l + V2)- The numbers in the figure indicate the values of the effective mass 
in the ground state at Fermi momentum for the calculations. The value of incompressibility is 
K = 380 MeV for the calculation represented by the long-dash-dotted line in the top panel and 
K = 210 MeV for other calculations in that panel. The solid and dotted lines in the bottom 
panel repeat the results from the top panel. The long-dash-double-dotted line in the bottom panel 
represents the results of a calculation where the momentum dependence of the MF at p > po is 
made to follow the dependence at p = po, cf. the text. 
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FIG. 7. Measured |^(| (filled symbols) and calculated (lines) ratios of out-of-plane to in-plane 
proton yields at midrapidity (0.35 < y/ybeam < 0.65), as a function of transverse momentum, in 
1000 MeV/nucleon 209 Bi + 209 Bi collisions. Both the measured and calculated ratios are computed 
from the lowest Fourier coefficients of the proton distribution, Rjy = (1— v-^/il+v^)- The top panel 
shows the ratios for collisions at b = 9.0 fm. The long-dashed, solid, and short-dashed lines represent 
calculations for the momentum-dependent MFs corresponding to K = 210 MeV and m* /m = 
0.65 (first of the sets in Table 0), m*/m = 0.70, and 0.79, respectively. The long-dash-dotted, 
dotted, and short-dash-dotted lines represent two calculations with momentum-independent MFs 
corresponding to K = 380 MeV, and K = 210 MeV , and one with no MF at all, respectively. 
The numbers in the panel indicate values of the effective mass for the calculations. The bottom 
panel shows ratios for collisions at impact parameters b = 11.0 fm (filled circles and dotted line), 
9.0 fm (filled squares and solid line), and 6.7 fm (filled diamonds and long-dash-double-dot line). 
The calculations for the bottom panel were carried out using the MF characterized by K = 210 MeV 
and m* /m = 0.70. The numbers in the bottom panel indicate the values of impact parameter. 
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FIG. 8. Time dependence of baryon density (stars in the top panel) at the system center and 
time dependence of midrapidity emission-rate for high-momentum (lines in the top panel) and all 
(bottom panel) protons, in Bi + Bi simulations at 400 MeV/nucleon and b = 8.7 fm. All results 
were obtained with the K = 210 MeV m* /m = 0.70 MF, except for the emission rate for the 
high-momentum protons represented by the thick dashed line, obtained with the K = 380 MeV 
momentum-independent MF. The thin vertical dashed line at t = 30 fm/c is drawn to guide the 
eye. 
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FIG. 9. Contour plots of baryon density in the 400 MeV/nucleon Bi + Bi reaction at b = 8.7 fm 
in the reaction plane (bottom panels) and in the plane through the system center othogonal to the 
beam axis (top panels). The simulation was done using the MF characterized by K = 210 MeV 
and m*/m = 0.70. The density contours are shown at intervals of 0.2 in p/po, starting from 
p/po = 0.1. Every third contour, i.e. for p/po = 0.1, 0.7, and 1.3, is represented by a dashed line. 
The remaining contours are represented by solid lines. Numbers in the figure indicate time in fm/c. 
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FIG. 10. Time dependence of — V2 for protons at p 1 - > 550 MeV/c from 400 MeV/nucleon 
Bi + Bi simulations employing either the m*/m = 0.70 K = 210 MeV MF (open circles) or the 
K = 380 MeV momentum-independent MF (filled circles). The top panel shows the average values 
for protons emitted in the vicinity of a given time, while the bottom panel shows the value for all 
protons emitted up to that time. 
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FIG. 11. Contour plots of the optical potential field, Eq. (|2|), in the 400 MeV/nucleon 
b = 8.7 fm Bi + Bi reaction, within the plane through the system center perpendicular to the 
beam axis, at t = 10 fm/c, for nucleons with p 1 - = 600 MeV/c directed to the right of the 
plot. The three panels, from left to right of the plot respectively, show results for three differ- 
ent MFs: the standard MF characterized by m* /m = 0.70 and K = 210 MeV, the same MF with 
momentum-dependence frozen above po> an d the momentum-independent MF characterized by 
K = 310 MeV. The Coulomb and isospin contributions to the potential are ignored in every case. 
The contours are indicated for the potential value of —1 MeV and for positive and negative mul- 
tiples of 10 MeV. The contours for the positive and negative values are represented, respectively, 
by the solid and dashed lines. The —1, ±10, and ±40 MeV contours are labeled. 
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FIG. 12. Proton Fourier coefficient t>4 = (cos 4</>) as a function of transverse momentum at 
midrapidity in the Bi + Bi system at 400 MeV/nucleon and b = 8.7 fm, from simulations utilizing 
the m* /m = 0.70 K = 210 MeV MF (solid line) and the K = 380 MeV momentum-independent 
MF (dash-dotted line). Numbers in the figure indicate the effective mass. 
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FIG. 13. Measured |4(| (filled squares) and calculated (lines) ratios of out-of-plane to in-plane 
proton yields at midrapidity (0.35 < y/ybeam < 0.65), as a function of transverse momentum, 
in 209 Bi + 209 Bi collisions at b = 8.6 fm and energies of 400 (top panel) and 700 MeV/nucleon 
(bottom panel). Both the measured and calculated ratios are computed from the lowest Fourier 
coefficients of the proton distribution, Rn = (1 — i>2)/(l + 1*2)- The calculations were carried out 
using in-medium cross sections and MFs characterized by the incompressibility of K = 210 MeV. 
The numbers in the figure indicate the values of the effective mass for the MFs; the solid lines rep- 
resent results obtained using the MF in Table | characterized by m* jra = 0.70 and the long-dashed 
and long-dash-double-dotted lines represent results obtained the first and second m* /m = 0.65 MF 
in the table, respectively. 
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FIG. 14. Optical potential in nuclear matter as a function of nucleon momentum, at different 
densities, from the DBHF calculations with Bonn-A interaction |43| and in our parametrizations. 
The short-dashed and short-dash-dotted lines represent the DBHF potentials obtained assum- 
ing, respectively, momentum-independent vector and scalar MFs [@,^3| and parametrized momen- 
tum-dependent MFs ||. The solid and long-dashed lines represent the optical potentials for our 
m* jm = 0.70 and first m* /m = 0.65 K = 210 MeV MF parametrizations, respectively. The num- 
bers in the figure indicate the values of density in units of the normal density po = 0.16 fm~ 3 . 
The thin horizontal dashed lines indicate the zero value for the potential. 
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FIG. 15. Optical potential in nuclear matter as a function of nucleon energy, at differ- 
ent densities, from the BBG calculations with Paris interaction (6|,^] and in our parametriza- 
tions. The short-dashed, dotted, and solid lines represent the BBG potential and the potential 
for our second m*/m = 0.65 and m*/m = 0.70 K = 210 MeV parametrizations, respectively. 
The long-dash-double-dotted lines represent the potential for our m* /m = 0.70 parametrization 
with the momentum dependence frozen above po = 0.16 fm -3 . The numbers in the figure indicate 
the values of density in units of fm~ 3 . The thin horizontal dashed lines indicate the zero value for 
the potential. 
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FIG. 16. Optical potential in nuclear matter as a function of nucleon energy, at different den- 
sities, from the variational calculations with AV14 + UVII interactions || and in our parametriza- 
tions. The short-dashed, dotted, and solid lines represent the AV14 + UVII potential and the po- 
tential for our second m* /m = 0.65 and m*/m = 0.70 K = 210 MeV parametrizations, respectively. 
The long-dash-double-dotted lines represent the potential for our m* /m = 0.70 parametrization 
with the momentum dependence frozen above po = 0.16 fm~ 3 . The numbers in the figure indicate 
the values of density in units of fm~ 3 . The thin horizontal dashed lines indicate the zero value for 
the potential. 
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FIG. 17. Optical potential in nuclear matter as a function of nucleon energy, at different 
densities, from the variational calculations with UV14 + TNI and UV14 + UVII interactions || 
and in our parametrizations. The short-dashed and long-dash-dotted lines represent the UV14 + 
TNI and UV14 + UVII potentials, respectively. The solid, long-dashed, and dotted lines represent, 
respectively, the potentials for our m* /m = 0.70, and first and second m* /m = 0.65 K = 210 MeV 
MF parametrizations. The numbers in the figure indicate the values of density in units of fm~ 3 . 
The thin horizontal dashed lines indicate the zero value for the potential. 
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